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We study the dynamics of the phase behavior of a polymer blend in the presence of 
shear flow. By adopting a two fluid picture and using a generalization of the concept of 
material derivative, we construct kinetic equations that describe the phase behavior of 
polymer blends in the presence of external flow. A phenomenological form for the shear 
modulus for the blend is proposed. The study indicates that a nonlinear dependence of 
the shear modulus of the blend on the volume fraction of one of the species is crucial for 
a shift in the stability line to be induced by shear flow. 



I. INTRODUCTION 



The dynamics of the phase behavior in polymer mixtures under external flow fields has 
aroused great interest over the last two decades Jl], ||, |3|, §|, |5], || [7|, |], |9], [L0| . The motivation 
for these studies is twofold. First, the effect of viscoelasticity on the phase behavior of 
polymer mixtures can be directly detected for some macroscopic flows. Secondly, many 
industrial processes, such as extrusion and painting processes, generate shear flow fields 
in polymer solutions and melts. In order gain further insight into non-equilibrium phase 
transitions and to optimize these industrial processes, one must understand the phase 
behavior of polymer mixtures in the presence of a flow field. Experimentally, a number of 
groups have reported that the phase behavior of polymer mixtures can be dramatically 
changed by macroscopic flow fields. In particular, for polymer-solvent mixtures in the 
presence of a shear flow, a greatly enhanced turbidity has been observed at temperatures 
much higher than the equilibrium critical temperature 0, • To study the mechanism of 
the observed phenomenon in polymer solutions, a number of theoretical efforts have been 
made |5], §, ||, pT0[| , and it is now understood that a nonlinear concentration dependence of 
the shear modulus is crucial for an upward shift of the phase separation temperature flO| . 



The temperature shift is proportional to the square of the shear strength in the regime of 
weak shear || |T0 . 



Although some experiments have been carried out for polymer blends under external 
flow Q, theoretically, one knows very little about the phase behavior of this system. 
Doi and Onuki first discussed the Langevin equations describing the dynamics of phase 
separation of a polymer blend 0. However, to our knowledge, the full consequences of 
the equations have not been explored. Moreover, the original approach employed in the 
ground-breaking work of Doi and Onuki prevents these kinetic equations from reducing 
to the polymer-solution case. Here, we present an expanded derivation, in which the 
difference between the Lagrangian and Eulerian descriptions has been taken into account. 
Using a phenomenological form for the shear modulus of a polymer blend, we can carry 
out a linear stability analysis for the model and find that the equilibrium spinodal line 
can be shifted in a complicated fashion by the shear flow. The purpose of this paper is 
to report on these studies and describe the relevant techniques in detail. 

The system consists of two kinds of polymers with different degrees of polymerization, 
Na and Nb- The volume fraction of polymer A at space-time point (r, t) is denoted by 
4>A(r,t) = 4>(r,t), and the volume fraction of polymer B is then 0s(r, t) = 1 — 0(r, t). 
We make the assumption that monomers of both species have the same specific volume, 
which can be expressed as || [7|, § 

Pa{t, t) pB{r,t) 

0(r,t) 1-0M) P ' U 

where pA and ps are the respective mass densities of polymers A and B and p is the total 
mass density, which is a constant. This assumption is consistent with the incompressibility 
of the system. In the two-fluid picture || |7|, || |9j , the two species are moving with different 
velocities, so that both bulk flow of the fluid and mutual diffusion of the two species, 



accompanied by chain deformation, take place simultaneously. Our task is to study the 
effect of flow and chain deformation on the dynamics of the phase behavior of the system. 

The paper is organized as follows. In Sec. II, we first discuss appropriate material 
derivatives and then derive the kinetic equations for the fluid velocity v and volume 
fraction (f). In Sec. Ill, we perform a linear stability analysis of the model, from which the 
effect of flow on the domain of linear stability can be obtained. Finally, our conclusions 
are briefly summarized in Sec. IV. 

II. THE MODEL 

A polymer blend is a viscoelastic system, sharing features of an elastic continuum and 
a viscous fluid. The system is elastic, but it only has a "faded" memory. The system is 
viscous, but it can bear deformation on some time scales. To derive the kinetic equations 
for such a system, one must call on familiar methods in studies of deformable media, as 
well as techniques for viscous fluids. 

A. Material Derivatives 

In continuum mechanics, one uses both Lagrangian and Eulerian coordinates to de- 



scribe the motion of a material element |Tl], [12| . The Lagrangian coordinates, which can 



be denoted by a = {aj}, i = 1,2, 3, are used to label the material elements or "particles" 
in a reference configuration (usually the undeformed state), while the Eulerian coordi- 
nates, denoted by r = {r^}, i = 1,2, 3, are the coordinates of the particles in the current 
configuration. The two coordinates are related through the following equations 

r = r(a,t) , (2) 
a = a(r, t) . (3) 

The velocity of the material element that is currently located at the point r is defined as 
the time rate of change of its position 

(4 ) 

where the subscript a is used to emphasize the fact that the derivative is to be evaluated 
for a particular material element whose Lagrangian coordinates are a. This is the usual 
material derivative. In the Eulerian description, the material derivative of any property 
pertaining to the particle labeled by the Lagrangian coordinate a is given by 

where v is the velocity of the particle a at position r, given in Eq. (f|). For small deforma- 
tions in which the Lagrangian strain tensor and the Eulerian strain tensor have the same 
form, the strain tensor e is the symmetric part of the displacement gradient tensor, 




1 V(Ar) + V(Ar) T , (6) 



2 



where Ar = r — a, and the dagger superscript stands for the transposition of tensors. The 
material derivative of the strain tensor is given by fTT|j , 



De 
Dt 



I( V « + W) . (7) 



Finally, the principle of conservation of mass leads directly to the following well-known 
formula [JTTJ, [TJ] 

-J d 3 rp(r,t)Q(r,t) = J d 3 rp(r, t)j^Q(r, t) , (8) 

where Q(r, t) is any physical quantity per unit mass and p(r, t) is the mass density of the 
material. 

We now generalize the above concepts to the two-fluid picture of polymer blends. As 
usual, we choose the undeformed state as the reference configuration, in which each ma- 
terial particle is identified by its Lagrangian coordinates a. Since in the present situation 
there are two kinds of material particles moving with different velocities in the system, 
one should distinguish the Lagrangian coordinates for the two species. We denote the 
Lagrangian coordinates of the particles of polymers A and B by and a# respectively. 
Then Eq. (^) can be generalized as 

(9) 
(10) 

The physical meaning of Eqs. (^j and (|10|) is as follows. At each space-time point in 
the current configuration, there are two velocities, iu(r, t) and t>s(r, t), which will be 
acquired by the material particles passing through this point, depending on the type of 
material particles. That is, particles of polymer A pass the point with velocity tu( r > t), 
while particles of polymer B pass the same point with velocity Vs(r, t). The fluid velocity 
(average velocity) of this point is given by 

v(r,t) = <j>(r,t)v A (r,t) + [l-(l>(r,t)]vB(r,t) . (11) 

Naturally, corresponding to Eqs. (^) and (|T0|) , we may introduce two kinds of material 
derivatives in the system 




D , d 



Dt )A dt 
D d 



or + fA-V, (12) 



+ v B -V. (13) 



If we focus on particles of polymer A, the material derivative is given by Eq. (0); similarly, 
the material derivative for species B is given by Eq. (|T3"[). The essential point is that any 
difference between va, vb, and the center of mass velocity is due to mutual diffusion. A 



constitutive relation is required to fix the diffusion flux or, equivalently, va — v B - We will 
return to this point later. 

It is easy to evaluate the two material derivatives for some basic physical quantities, 
such as the volume fraction and strain tensor e. First, from the continuity equations 
for both species 



dp A (r, t) 
dt 

dpBjryt) 
dt 



+ V • p A (r, t)v A {r, t) = 0, 
+ V- p B (r,t)v B (r,t) = 0, 



we can obtain the expressions for the material derivatives of the volume fraction 

~D<j>(r,ty 

J A 

= [l-0(r,t)]V-Mr,t) 



Dt 

D<p(r,t) 



(r,t)V • v A (r,t) 



Dt 



(14) 
(15) 

(16) 
(17) 



B 



where Eq. ([[]) has been used. Note that augmented by Eq. ([!]), the continuity equations 
lead directly to the incompressibility condition 

V • v = V • <fw A + V • (1 - (p)v B = . (18) 

Next, for small deformations, the strain tensor of the polymer blend is still given by 
Eq. (OF), but the material derivative of e is generalized to the following equations, 



f- 

\Dt 



2 V~ 



A B . 



(19) 



Taking material derivatives defined in Eqs. flT2p and (|T3|) on both sides of Eq. (|5p, and 
noticing that for small deformation, V = d/dr ~ 8/ da, Eq. (|i~9D is obtained. As will 

and ([TOD are useful in the evaluation of the 



be seen in Subsection C, Eqs. (|l6j), (hj 
dissipation rate of the total free energy of the system. 

Finally, it is follows from the principle of mass conservation that Eq. (f|) still holds. 
But in the present situation, p is the total mass density p = p A + Pb, which is a constant 
(see Eq. (|l|)), and v is the fluid velocity given in Eq. (|TT1) . Noticing the fact that the 
system is incompressible (V • v = 0), for a polymer blend Eq. (§) can be written as 



dQ 



(20) 



Here, a boundary term has been ignored. Furthermore, since the masses of both species 
are also conserved individually, one has the following equations 



d 
dt 



d 3 rpiQi(r,t) 



d 3 rp r , 



(- 

\Dt 



Qi(r,t) 



A,B 



(21) 



where piQi is any physical quantity contributed by the z-species. Note that only two of 
the three formulae for the time derivative of the volume integration are independent. In 
fact, summing up the two equations given in Eq. fl2l), one obtains Eq. (||). 



B. Total Free Energy 

We take into account three kinds of contributions to the total free energy of the system: 
The kinetic energy K of moving particles, the mixing free energy F m of the two species, 
and the elastic free energy F e of polymers due to chain deformation ||. Thus, the total 
free energy of the system can be written as 

F t = K + F m (0) + F e ((f>,e) . (22) 

Here, we have assumed that the mixing free energy is a function of <fi only, while the 
elastic free energy depends on both and e. The kinetic energy of the two kinds of 
moving particles can be expressed as 

K = jd 3 r (±p A v 2 A + l -p B vl) . (23) 

The mixing free energy can be written as 

F m {4>) = ( d 3 rf m (<P) , (24) 



where f m {4>) is the mixing free energy density, which, for example, can be chosen to be 
the Flory-Huggins form. Our derivation for the kinetic equations is independent of the 
precise form of f m (<f>). In a similar way, the elastic free energy can be expressed as 

F e (0,e) = J d 3 rf e (<f>,e) . (25) 

In the theory of linear elasticity, the elastic energy density due to chain deformation can 



be phenomenologically expressed as P, WA 



/ e (0,e) = G(cf>)r e (e), (26) 

where /* (e) = e : e, and the coefficient G is the shear modulus of the system, which, 
in general, depends on concentration only. Here, the notation (:) stands for the scalar 
product of second rank tensors. 

Since in a polymer blend both species contribute to the elastic energy, we propose the 
following intuitive form for the shear modulus of the blend 



G = Gg> 



n (o) r (o) 



A(0) , (27) 



where £?f^ is the shear modulus of the i-species before mixing ("bare" shear modulus), 
and A(4>) is an interpolating function describing the effect of blending. The condition 
that G((f> = 0) = G { b } and G(<p = 1) = Gf requires that A(0) = and A(l) = 1. Then 
the simplest form of the interpolating function would be A(0) = 0, i.e., an ideal mixture 
approximation. More generally we suppose 



A(0) = <p +A(0) 



(28) 



with A(0) = A(l) = 0. For use in Eq. (|26|) , it is convenient to cast Eq. 

G{<j>) = <PG A (4>) + (1-0) G B (<f>) , 



into the form 
(29) 



where Gi, i = A,B, are the "renormalized" shear moduli of species A and B individually, 
which can be expressed as 



Ga{<P) 

g b {<I>) 



A 



1 + 



A(<j>) 



A(0) 
1-6 



(30) 
(31) 



As we will see below, Eq. (p9~D is a reasonable approximation leading to a sensible result 
for the network velocity. 

Finally, since the polymer blend is viscoelastic (rather than a purely elastic system), 
one expects that Eq. (^) describes the physics within the time scales less than the relax- 
ation time of the shear stress (usually referred to as the "terminal relaxation time" ) |T7[ . 

C. Dissipation Rate of Total Free Energy 

With the results presented in the previous two subsections, we are ready to discuss 
the dissipation rate of the total free energy of the system. First, using the formula given 
in Eq. (|2"TD , the evaluation of the time derivative of the kinetic energy is straightforward, 
and the result can be expressed as 



d 4 = f d s r L AVA .(^) 



dt 



■Dva 



\ Dt ) A 



+ Pbv b 



(Dv 



B 



V Dt 



B 



(32) 



Next, since, in general, the mixing free energy f m {4>) cannot be simply divided into contri- 
butions by particles A and B, it is convenient to use Eq. (p0|) to calculate the dissipation 
rate of the mixing free energy. Indeed, it is easy to obtain 



3 df m (<j>) 



d 6 r 



dt 



3 df m (4>) d(p 



d 6 r 



dt 



(33) 



Making use of the continuity equation (|lq) and integrating by parts, we have 



dF m 
dt 



J d 3 r v a ■ 0V 



dfm((f>) 



J d 3 r v a ■ V7r r , 



(34) 



where 7r m is the osmotic pressure associated with the mixing free energy, given by 



7T, 



(35) 



Note that if the continuity equation ( |T7| ) were used, a different expression for dF m /dt 
would be obtained, but it will give the same final kinetic equations when the condition 
V • v — is taken into account. 



Finally, we discuss the time derivative of the elastic free energy. In view of Eq. ((29]), 
we can write 



fe 



/ eA (0,e) + f eB ((f),e) 



(36) 



where f e i((fi, e) = 0iCj/*(e), for i = A, B, are the elastic free energy densities contributed 
by the two components individually. Using Eq. (^1|), we have 



dt 



d 3 r \ (j) 



jj;4> VeA(0,e) 



+ (1 



*(1 

DV 



feB{<f>,e) 



■ (37) 



Using the chain rule, the material derivatives of the free energies f eA (<fi,e) and f eB (4>,e) 
can be expressed in terms of the material derivatives of and e. Making use of Eqs. (|16"D, 
(fnp, and fllTf), we obtain 



5/, 



d d r <j -vr eA V ■ w A - vr eB V ■ u B + (Vv A ) : + (Vv fl ) : 



B 



where 7r ei , with % = A, B, are the "elastic osmotic pressures", given by 



7T fi 



^-lVei(0,e) ,i = A,B 



(3? 



(39) 



The stress tensor acting on the network can be defined as [[L5 



T 



In view of Eqs. (|S^) and (PB"D, it is easy to check that 

dfeA , Ca 



<9e 

dfeB 

dt 



G 



T 



B 



G 
] ~G 



T . 



(40) 

(41) 
(42) 



Substituting Eqs. ((41]) and (|42| ) into Eq. (^), the dissipation rate of the elastic free energy 
can be written as 



dFe 
dt 



d 3 r \ Va 



Vvr eyl - V • 0^r r 



+ v B 



, (43) 



after an integration by parts. 

Combining Eqs. fl32|), (|34"D , and (|43"D, we obtain the dissipation rate of the total free 
energy 

dF t r 3 r /£>u A \ /£>^b\ 

~dt = l dr \ pAVA \-m) A + pBVB \-Dt) B 



+ va ■ 
+ v B ■ 



V(7T m + 7l eA ) - V • 0^T r 



V7T eB - V ■ (1 



<G 



-T 



(44) 



It should be understood that the partial derivative d/d(p is carried out at fixed e and d/de 
is carried out at fixed 0. 



D. Network Velocity 

In this subsection, we discuss the network velocity (or tube velocity in the reptation 
picture). Substituting Eqs. ( f4"T| ) and ([4^) into Eq. fl38]), the dissipation rate for the elastic 
free energy can be expressed as 



dFe 
dt 



J d 3 r {-n eA V ■ v A - 7r eB V • v B 

G- 1 [(f>G A Vv A + (1 - <f>)G B Vv B ] : r} . (45) 



The last term in the above equation describes the time rate of change of the elastic energy 
purely due to the change of the strain tensor, so that the coefficient of the stress tensor r 
can be identified as the gradient of the network velocity (or tube velocity in the reptation 
picture), 

Vvt = G- 1 [<j)G A Vv A + (1 - <P)G B Wv B ] . (46) 

As we now show, this expression is in agreement with a previously obtained result using 
a microscopic approach. 

The network velocity v t has previously been estimated from molecular theory [7|, [16[ 
for a uniform system (0 is constant) with the result 

(ava + (bv b . . 

Here Q, i — A,B, are effective friction coefficients given by 

Ci = ^i^Co, i = A,B, (48) 

where vi is the number of chains of species % per unit volume, Nf is the average interval 
between two successive entanglement points along one chain, and (o is the phenomenologi- 
cal friction coefficient between the two species. On the other hand, when is independent 
of space, Eq. ( fi6|) can be solved with the result 



v t 



G- 1 [<PG a v a + (1 - <P)G B v B ] , (49) 



where an integrating constant has been determined as zero from the condition that v% = 
v A , when = 1. We may suppose that G A)B measure the densities of entanglement points, 
i.e., Gi oc T0;/iV7 0. Then Eq. (0) is recovered. 

E. Kinetic Equations 

First we derive the equations for the two velocity fields v A and v B . These equations 
can be obtained by means of Rayleigh's variational principle || |7J. Following Doi and 
Onuki p, [/J, one can define a Rayleighian functional, 



where F t is the total free energy of the system and W is the dissipation function due to 
relative motion of the two polymers, which one assumes can be written as 



W 



d 3 r c(r)C (v A - v B f 



(51) 



Here c is the monomer concentration of species A defined via = v m c with v m the 
monomer volume, and ( is the friction constant, which, in general, is a function of <f) ||. 
The Rayleighian defined in Eq. fl5"U| ) may be understood as the total energy dissipation 
rate of the system. The variational principle states that the velocities v A and v B are 
determined by the condition that they minimize the Rayleighian || [7|. 

Substituting Eqs. (|44l ) and (^l|) into Eq. fl50f) , we have the Rayleighian for the polymer 
blend 



R 



PaVa 



(Dv A \ 



+ va ■ 



+ v B ■ 



V(vr m + 7T eA ) - V 
VTTeU - V • (1 - q 



+ 



^c(r)C {v A -v B ) 2 } 



+ Pbv b ■ 
Gb 



(Dv B \ 



-T 



(52) 



Since v A and v B are not independent variables due to the incompressibility condition ([18]) , 
the functional R must be minimized under this condition. The conditional minimization 
of the functional R with respect to v A and vb leads to the following equations, 

G A 



(Dv A \ 



vb 



PB 



(Dvb\ 
\ Dt A 



-c({v A 
cC,{v A - v B ) 



l-0)(Vp) 



V(7T„ 



T^eA) 



G 



■T 



Vn eB + V • (1 - 



G 



■T 



(53) 
(54) 



where p is the Lagrange multiplier imposing the incompressibility condition. Eqs. (|53|) 
and fl5"4|) describe the motion of polymers A and B in the system. 

Solving Eqs. (|5^) and (0) for v A and then substituting the resulting expression into 
Eq. (0), we can obtain the diffusion equation for 0. Clearly this cannot be done exactly 
and some approximation must be applied J7|. Since the velocities relax much faster than 
<f), we can ignore the inertia terms in Eqs. (^) and (|54]) to obtain an explicit expression 
for v A . After eliminating the p terms, we can express v A as 



v|£-a[(Ve) 



: t + V 



T 



P 



T 



(55) 



where / = f m + f e and Eq. (|TT|) has been used. In Eq. (0), we have introduced two 
parameters a(0) and (3(4>) for convenience, which are given by 



P = 



G 1 — G_b] 
G A 



V(l 



(56) 
(57) 



Substituting Eq. ([55]) into Eq. (0), we obtain 

1 - 



d( t> 



c 



- { V%- - a [(Ve) : r + V • r] - /? • r) , (5? 



where the incompressibility condition has been used. This is the diffusion equation for 0. 

It is convenient to describe the motion of the system using the fluid velocity v and 
relative velocity u = va — vb- From Eqs. (^) and fl54|) , it is easy to check that the kinetic 
equations for v and u can be written respectively as 



-Vp - V7r + V • r 



(59) 



1 - 



<X M -^i V ^T -«[(Ve) : r + V • r] -/3-rl (60) 



where 7r = 7r m + 7r ej4 + 7r eB is the total osmotic pressure, and R v and R u are couplings 
between v and u that can be expressed as 



R v 
R u 



)u + 



)u ■ Vu 



u ■ V(u — du) + (v — <bu) ■ Vm 



(61) 
(62) 



Eq. (|59D is the generalized Navier-Stokes equation for the fluid velocity. R v is the cor- 
rection due to coupling between bulk flow and relative motion. Eq. (|60|) describes the 
relative motion between the two species. 

The constitutive equation cannot be derived from the above formalism, because so far 
no ordinary viscosity effects have been included in the Rayleighian. We assume that the 
time evolution of the stress tensor is described by the upper convected Maxwell equa- 
tion 0, H, [18 



A 



da 



Vcr - cr • Vv t - (Vff) 1 " • cr 



a 



G(<P)6 , 



(63) 



where A is the relaxation time, v t is the network velocity, and the stress tensor a is 
related to r through a = G5 + r. The network velocity v t has been used here because 
the stress acts on the polymer network |7| [19]. Since in polymer blends the fluctuations 
of concentration are small, and Eq. ([R]) cannot be solved explicitly, we will use Eq. ( ^9|) 
as an approximate expression for the network velocity Vt in the linear analysis. Using the 
expressions for va and Vb, we have finally 



v - v m a- 



a [(Ve) : r + V • r] - j3 • r 



(64) 



Eqs. ([58]), (|59]), (|5UD, and ( [53] ) describe the dynamics of phase behavior of polymer blend 
under shear flow. Since u appears only in Eqs. ( |59"P and (|B"UP, and in most cases R v is 
negligibly small 0, one can ignore the R v term in Eq. (^9]), and then the model reduces to 
three equations. In the case where the S-component represents small molecules, Gb ~ 0, 
a ~ 1/0, (3 ~ 0, the model for a polymer solution is recovered @, TTJ . 



Note that although the structure of our set of kinetic equations is similar to that given 
in Ref. J7[, some important differences exist. First, the free energy density / in Eq. ( |58"D 
is the total free energy density, i.e., f = f m + f e , while in Eq. (4.17) of Ref. [0], / is the 
mixing free energy only. That is, we have taken into account the fact that, in general, the 
elastic free energy can be dependent (through the shear modulus G(<ftj). Second, two 
more terms (a(Ve) : r) and {[3 ■ r) are generated in Eq. d58|). The a(Ve) : r) term and 
the inclusion of f e in the total free energy density / ensure the that kinetic equations to 
reduce appropriately to the solution case. 



III. LINEAR ANALYSIS 



Eqs. (0), (0), and fl63|) are nonlinear equations, so that a complete analysis is difficult. 
However, to study the effect of viscoelasticity on the phase boundary under shear, it is 
sufficient to carry out a linear analysis similar to that developed in Refs. |5|, ||. In this 
formalism, we first iteratively solve Eq. (F3) for a to the "second-order fluid" level from 



which the stress tensor a can be expressed in terms of and v \TE]. Then we substitute 



this constitutive relation for a(v,(f>) into Eqs. ( p8|) and (|59|) . Setting 

= 0o + 0i , (65) 
v = v + v i , (66) 

where 0o and v o are the overall average volume fraction and the fluid velocity, and 0i and 
Vi are small deviations, we can solve the equations to linear order in X and V\. In the 
case of shear flow, vq = Sye x , where S is the shear rate and e x is the unit vector in the 
x-direction. The expansion parameter is essentially the combination SX, which limits the 
approximations to the regime of weak shear. 

A. Iterations for o and e 



We first solve the Maxwell equation iteratively |18[]. Since our aim is to see how 



viscoelasticity changes the stability line for phase separation, it is sufficient to obtain 
a solution for a in the long wavelength limit. It is easy to check by "power counting" 
that 0i and Vt>i are the leading order terms, and then v — v t ~ V0 + V 2 v are higher 
order terms. Therefore, we can replace Vt by v in the Maxwell equation for the present 
purpose. Furthermore, since v and a relax much faster than 0, we may ignore the inertia 
terms p, [| and obtain 

A [-a ■ Vv - (Vv^ -a]+a = G(0)5 . (67) 

Directly iterating to the "second-order fluid" (i.e., to 0{\Vv | 2 ), we obtain 

a = G5+ 2r]D + 2r]\ (d ■ X7v + W ■ D) , (68) 

where i] = XG is the viscosity and D is the gradient tensor given by 2D = Vv + X?v' . 
Note that, at this order, only leading order non- Newtonian terms are included. Finally, 



for later convenience, we write down the expression for the strain tensor e 

e = XD+ X 2 (D ■ Vv + W ■ D) . (69) 



Now we have expressed a in terms of 4> and v. Our next task is to eliminate a and e using 
Eqs. fl5"B|) and (^) in Eqs. (|5*%D and (|59|) and carry out a linear analysis for the remaining 
equations. 

B. The Navier-Stokes Equation 

Now we discuss the linear analysis of the Navier-Stokes equation. In a similar spirit to 
the treatment of the Maxwell equation, we ignore the inertia terms since fluid velocities 
relax much faster than [||, ||. Then the equation for v becomes simply 

Vp + Vtt - V • r = . (70) 

We next find expressions for Vtt and V ■ r linear in (f>i and v±. For later convenience, 
we introduce two parameters g and £ through G(<p) = /csTg(0) and rj(4>) = /cbT£(0). It 
follows from the definition that it = 7r(</>, e) = 7r(0, v ), and therefore we have 

V.. (|) o V, 1 + (V.).(f) o , (H) 

where the subscript "0" indicates that the derivatives are evaluated at 4>o and vq. In view 
of Eq. (|68D , we can express V ■ r in a similar way as 

V • t = k B TV ■ [2£D + 2£A (p ■ Vv + Vv^ ■ £>)] . (72) 

Substituting Eqs. (|65D and (|66|) into Eq. ([72]) and keeping terms only linear in 0i and v\, 
we have 



V-r = knT 



e V 2 t;i + (e^^i + eydM + (9(^ 2 , 5^)1 , (73) 



where with % = x,y,z are unit vectors. Here the subscript "0" means that the values 
are evaluated at <fi = 0o? an d primes indicate ^-differentiation. To obtain the shift in the 
spinodal (stability line) to leading order in the shear rate (O(SX) 2 ), the 0(S 2 , Sv) terms 
can be ignored at this stage. 

Substituting Eqs. (ffT|) and (|73| ) into Eq. (|70|) and eliminating p via V • v — 0, we can 
solve f to order 0(5*) in Fourier space with the result 

i€o k y ( n i2 



v lx (k) = -iS^ {2k' x - 1) Mk) , (74) 

vi y (k) = -^1^(2^-1)^), (75) 

v lz (k) = -2iS&!^Mk), (76) 

where ki = hjk. As we see, to this order, the solution for v is independent of the precise 
form of 7r. Note that the expressions given in Eqs. are the same as those for 



polymer solutions || [10| . That is, to leading order, the linear <f> — v relations are the same 
for both polymer solutions and blends. 



C. Shift in the Stability Line 



We now discuss the linearization of the diffusion equation (0), from which the effect 
of viscoelasticity on the phase boundary can be studied. We choose the Flory-Huggins 
form for the mixing free energy p 



fm((f>) 



k B T 



N A v N B 



ln(l - 4>) + X <f>(l ~ 0) 



(77) 



where x is the Flory-Huggins interaction parameter, and use Eq. (F2^) for the elastic free 
energy. Substituting Eqs. ( |S5] ) and ( pED into Eq. fl5"8"| ) and making use of Eqs. (|7^)-([7oT), 
we can obtain the following linearized diffusion equation in Fourier space 



- - Sk x — 
dt x dk„ 



kElT a n 

-V m — <Po{l 



?k 2 



y + y 1 (k x ,k y ) Mk) ■ (78) 



Here ^ is an isotropic constant given by 

* O (0o) ' 1 



2 X + k(0 o ) (SXoY 



N A <f) N B (l-<fa) 

where Aq = £,o/go is the average relaxation time and the coefficient k is given by 



(79) 



k(0o) 



-a"- 9 l 
2 % go 



"Ofi'o 



JO) 



Note that k depends only on the shear modulus of the blend. The second term on the right- 
hand-side of Eq. ( |78D produces an anisotropic modification of the scattering function H 
and is given by 



= -2S\ (A + B )k x k v + (S\) 2 [4(Co-2B Q )(k x k y ) 2 
- 2(A + B - Do)k 2 x + (Co - Bo)k 2 } 

where A , B , Co, and D are the constants 



Bo 

Co 
Do 



dodo + <Po 1 9ao + (1 - <f>o) 1 9bo , 
ao9o-z- , 

so 
aog'o- 



(81) 

(82) 
(83) 

(84) 
(85) 



The zero-shear spinodal line (linear stability line) in the X ~ parameter space is deter- 
mined by the first three terms of [|17| . The last term is the modification of the stability 
line arising from the dynamics of viscoelasticity; this shift is 0(S 2 ), where S is the shear 
rate. 



The direction of the stability line shift is dependent on the sign of k, which in turn 
depends on the details of the interpolating function A(</>). When k > 0, the effective 
value of x is reduced, and stability line in the x ~ 4> plane is shifted to larger values of 
X (lower temperatures). On the other hand, when k < 0, the stability line is shifted in 
the opposite manner. Since, in general, the shear moduli for the individual species A 
and B involve material parameters, the details of the shift in the stability line cannot be 
determined explicitly. In highly symmetric situations, the magnitude of any shift might 
be small, owing to the small value of the coefficient [G^ — G^]. However, the shift in the 
linear stability should, according to the present analysis, correlate with measurements of 
the interpolating function for the blend, A (</>). If A(0) = <ft, k = and there is no shift 
in stability. 



IV. CONCLUSION 



We have applied the general two-fluid approach of Doi and Onuki ]7) to establish 
model equations for the study of the phase behavior, or more precisely, linear stability, 
of a polymer blend in the presence of shear flow. A phenomenological form of the shear 
modulus for a polymer blend is used, which can be experimentally determined. The 
modifications used here allow the kinetic equations to reduce to the solution case in the 
limit in which one of the species becomes a small (solvent) molecule. Linear analysis of the 
model indicates that the equilibrium stability line is shifted by the effect of viscoelasticity 
when the dependence of the shear modulus of the blend on the volume fraction of one of 
the species is nonlinear. The direction of the shift is dependent on the material properties 
of the species and on the range of volume fraction, so that the nature of the temperature 
shift is more complicated than that in polymer solutions |8], 1Q| . Physically, this feature 
of polymer blends can be attributed to the fact that both species have contributions to 
the elastic free energy. 
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